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Abstract— Image segmentation is the first step of image 
processing and image analysis. Texture segmentation is a chal- 
lenging task in image segmentation applications. Neutrosophy 
has a natural ability to handle the indeterminate information. In 
this work, we investigate the texture image segmentation based 
on Gabor filters (GFs) and neutrosophic graph cut (NGC). We 
proposed an image segmentation approach, which applies GFs 
to gray-level images to extract image features matrix, and it 
segments them into regions. First, color images are transformed 
to gray level images as input images. Then, input parameters of 
GFs are adjusted, and GFs are performed on the input images 
to extract features. The NGC is employed for classification of 
input images. Finally, experiments are conducted on various 
natural images to evaluate the approach. Experimental results 
show that the proposed approach achieves desired performance 
of texture segmentation. However, it cannot segment the texture- 
free images as well as texture images. In future works, we will 
try to segment both texture images and texture-free images at 
the same time. 
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I. INTRODUCTION 


The goal of the image segmentation considered as dividing 
an input image into several sub-images. To do this, some 
predefined criterions should be accepted where the sub- 
images are disjoint, homogeneous and meaningful. Many 
researchers have been working on image segmentation, and 
a variety of works have been published [1]. Texture seg- 
mentation is a challenging and important task in image 
segmentation applications. For the texture feature extraction, 
there exist various texture descriptors such as gray level co- 
occurence matrix, Gabor filters (GFs) and Markov random 
fields [2]. GFs have been successfully applied to many image 
applications such as texture segmentation [3-5], document 
analysis [6], edge detection [7], retina identification [8], 
fingerprint processing [9], image coding [10], and image 
representation [11]. 

Graph based methods are important for image segmen- 
tation applications [12]. A graph G can be defined as 
G = (V, E) where V and E are a set of vertices and edges, 
respectively. On an image, vertices can be either pixels or 
regions, and edges connect the neighbouring vertices [13]. 
A weight, which is a non-negative measure of dissimilarity, 


is associated with each edge based on some property of 
the pixels that it connects. Normalized-cut, Minimal-cut and 
Mean-Cut can be seen in some of popular graph cut methods 
14-16]. Neutrosophic set (NS) was proposed as an extension 
of fuzzy set [17]. In NS theory, a member of a set has 
a degree of truth, a falsity degree and an indeterminacy 
degree [18]. Therefore, it has an ability to deal with the 
indeterminacy information, and has attracted much attention 
almost in all engineering communities and subsequently a 
great number of works have been studied [19], [20]. 

In this paper, we propose a texture image segmentation 
algorithm which uses GFs, NS and graph cut method. Firstly, 
various GFs are considered which constructs a texture feature 
volume. Constructed volume is then projected into 2d by 
using principal component analysis (PCA). The projected 
feature image is transformed into NS domain and indeter- 
minacy degree is calculated accordingly. Then an indeter- 
minacy filter is constructed using the indeterminacy value 
on feature image which is defined by combining the spatial 
information and global texture feature information. The 
indeterminacy filter is employed to remove the indeterminacy 
in the intensity and spatial information respectively. A graph 
is defined on the image and the weight for each pixel is 
represented using the value after indeterminacy filtering, and 
the energy function is also redefined using the neutrosophic 
value. Finally, the segmentation results are obtained using a 
maximum-flow algorithm on the graph. 

The organisation of the paper is as following; in Section 
2 we give the methodology of the proposed method. In 
Section 3, some experimental works are conducted. Finally, 
we concluded the paper in Section 4. 


II. METHODOLOGY 


Briefly speaking, the proposed approach composed of 3 
main building blocks. Extracting texture features from input 
image by using GFs is the first building block. This block 
constructs a texture feature volume according to the number 
of considered GFs. The second block employes the NS and 
calculates the indeterminacy degree accordingly. Then an 
indeterminacy filter is constructed using the indeterminacy 
value on feature image which is defined by combining the 


spatial information and global texture feature information. 
This filter aims to remove the indeterminacy in the texture 
feature and spatial information respectively. Finally, a graph 
is defined on the image and the weight for each pixel is 
represented using the value after indeterminacy filtering, and 
the energy function is also redefined using the neutrosophic 
value. Finally, the segmentation results are obtained using a 
maximum-flow algorithm on the graph. The flowchart of the 
approach is illustrated in Fig. 1. 
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Fig. 1. The flowchart of the approach. 


The GFs theory was described in Jain and Farrokhnia [21]. 
GFs should be in array form and it is designed with different 
orientations and frequencies. The step of orientation angle 
can be 30, 45, etc. degrees and orientations are shown as 
below: 

Orients = [0,a,2 x a,---,180— al (1) 


where a is the step of orientation angle. Wavelengths have 
values in increasing powers of two start from minimum 
wavelength to maximum wavelength. 
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where Wmin is the minimum value of wavelength, Wmaz iS 
the maximum value of wavelength. 1m,;ows and imecois are 
the size of rows and columns of the input image. 
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Waves = [2°,2',2?,---, Be?)) X Wmin (5) 


The set of orientations and wavelengths are taken from [21]. 
Working GFs mean working with Gabor magnitude response 
of each filter. The size of the magnitude response is: 


Orientsien X WaveSien (6) 


where Orientsien is the length of array of orientation 
Orients in Equation (1), and Wavesje,, is the length of 
array of wavelength Waves in Equation (5). 

When using Gabor magnitude responses as features, some 
post-processing is required. Feature smoothing and normal- 
ization is applied as post-processing. In smoothing process, 
the Gaussian process and normalization is applied in order 
to get zero mean and unit variance. Smoothing is important 
because some local variations in each Gabor response are 
not suitable for the segmentation and they can be soften 
by using a low-pass filtering. In addition, spatial x and 
y axis properties of each pixel are added to the features 
matrix. Spatial location information improves the segmen- 
tation process to prefer groupings that are close together 
spatially. After the post-processing, we apply the PCA in 
order to reduce dimension size of the feature vector. PCA, is 
a procedure which uses an orthogonal transformation [22]. 
The transformation is carried out with the first principle 
component. 

After PCA, the reduced feature matrix is fed into the 
NS domain. A neutrosophic image [yg is interpreted using 
three subsets 7,, J, and F,. Based on the texture feature 
value and local spatial information, the true membership 
and indeterminacy membership are used to describe the 
indeterminacy among local neighbourhood as: 


g(a, y) — Gmin 


Imax — 


Ts(2,y) = (7) 


Gmin 
Gd(z,y) — Gdinin 


where g(x,y) and Gd(z, y) are the gray scale texture feature 
value and gradient values at the pixel of (a, y) on the image, 
respectively. The algorithm also computes the neutrosophic 
membership values based on the global texture feature value 
distribution which considers the indeterminacy on texture 
feature value between different groups. To this end, NCM 
is used to obtain the indeterminacy values between different 
groups on texture feature value to be segmented [23]. Before 
calculation of the global texture feature value based NS 
memberships, the obtained local Iyg is filtered by inde- 
terminate filter based on J,. After, constructing Iyg, an 
energy function for graph model is determined and the Ing 
is partitioned using maximum flow algorithm [24]. 

The algorithm of the proposed method is as following: 

Step 1: Compute the GFs of input texture; 

Step 2: Smoothing and normalization Gabor responses; 

Step 3: Dimensionality reduction with PCA; 

Step 4: Compute the local neutrosophic value T; and Is; 

Step 5: Filter J’, using indeterminate filter based on J,; 

Step 6: Employ NCM algorithm on the filtered T, subset 
to calculate 7;, and I,,; 

Step 7: Filter T;, using indeterminate filter based on I,,; 

Step 8: Define the energy function based on the T’, value; 


(8) 


Step 9: Partition the image using the maximum flow al- 
gorithm. 


III. EXPERIMENTS 


Experiments are conducted on various natural images 
from Berkeley Segmentation Dataset (BSD) to evaluate the 
approach. Four images are selected randomly to show the 
performance of the method. These images are “Cheetah”, 
“Lady”, “Plane” and “Zebra”, respectively. As shown in Fig. 
2, each of figure which includes these images has four types 
of images. First image is original gray-level and second 
image is after Gabor filtered and PCA image. Third image 
is segmentation result of NGC of the first. Fourth image is 
border result of NGC of the first. 


Fig. 2. Top-left: Cheetah image, top-right: result image after GF and PCA, 
bottom-left: segmentation of NGC, bottom-right: border of segmentation. 


GFs are applied to the input images to extract texture 
features. To use GFs as features, the orientations and the 
wavelengths should be selected and calculated. We consider 
the “Cheetah” image in Fig. 2 to examine the process of 
feature extraction. As seen in Fig. 2, the input image is 
a gray-level image. Its size is 321x481 pixels. Minimum 
wavelength is mentioned at upper section and it is 2.8284. 
Maximum wavelength is calculated via size of the input 
image. Maximum wavelength is 321? + 481? = 578.2750. 
The wavelengths that send to Gabor magnitude as input 
parameter, start in increasing powers of two from minimum 
wavelength up to maximum wavelength. For the “Cheetah” 
image, the wavelengths have six values: {2.8284, 5.6569, 
11.3137, 22.6274, 45.2548 and 90.5097}. 

Orientation angle called theta begins with 0 degree and 
continues with 45 degrees increase thus it has four values 
for the “Cheetah” image: {0, 45, 90 and 135}. 

The related illustrations are given in Figs 2-5. As can 
be seen in these figures, the proposed method produced 
reasonable segmentation results. The objects in the input 
images are well segmented from the background. However, 
the boundaries of the object are rough and need further post- 
processing. Especially, for woman image, the head of the 
woman is not correctly segmented. In addition, several over- 
segmentations are clear. 


Fig. 3. Top-left: Lady image, top-right: result image after GF and PCA, 
bottom-left: segmentation of NGC, bottom-right: border of segmentation. 


Fig. 4. Top-left: Plane image, top-right: result image after GF and PCA, 
bottom-left: segmentation of NGC, bottom-right: border of segmentation. 


IV. CONCLUSIONS 


This paper presents a texture image segmentation using 
GFs and NGC method. Before the image is described in 
NS domain, GFs are applied on the image to extract texture 
features. The experimental results show that the NGC method 
can segment the images properly and the proposed approach 
achieves desired performance of texture segmentation. The 
disadvantage of the proposed method is rough segmentation. 
To overcome this disadvantage, in the future works, a bound- 
ary preserving filtering method will be adopted. 
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